Revised Fluxes for Gliese 229B 



S.K. Leggett 

Joint Astronomy Centre, 660 N. A 'ohoku Place, Hilo, HI 96720 
skWjach. Hawaii, edu 

D.W. Toomey 

NASA IRTF, P.O. Box 4729, Hilo, HI 96720 
toomey@irtf.ifa.hawaii.edu 

T.R. Geballe 

Gemini North Observatory, 670 N. A 'ohoku Place, Hilo, HI 96720 
tgeballe@gemini. edu 

and 
R.H. Brown 

Univ. of Arizona Lunar & Planetary Lab., Tucson AZ 85721-0092 
rhb @lpl. arizona. edu 

Accepted March 25th 1999 for ApJ Letters 



ABSTRACT 

We have used the coronographic instrument CoCo with the IRTF's facility camera 
NSFCAM to obtain improved photometry at JHKL' of the giant planet/brown dwarf 
Gliese 229B. We have recalibrated the published spectra for this object, and re-calculated 
its luminosity. Our L' value, and our flux calibration of the spectra at JHK, are signif- 
icantly different from those previously published. Our results show good agreement at 
all bands except H with evolutionary models by Burrows et al. which include grain con- 
densation. The model comparison implies that Gliese 229B is likely to be a 0.5 Gyr-old 
25 Mj object with T eff -900 K. 

Subject headings: stars: individual (Gliese 229B) — stars: low-mass, brown dwarfs 
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Introduction 



Nakajima et al. 1995 discovered Gliese 229B dur- 
ing a coronographic search for low-mass compan- 
ions to young M-dwarfs. Despite continued searches, 
Gliese 229B to this date is the only substellar object 
that is resolvable and cool enough to show planet-like 
absorption features due to methane. A good under- 
standing of this object is key to our understanding 
of brown dwarfs and extra-solar planets in general. 
Published work on this unique object have been pri- 



Geballe et al 



1996 


pppenheimer et al. 1998) or theoretical (e.g. 


Allard et al. 1996, Marley et al. 1996|, fTsuji et al. 


1996 


Burrows et al. 1997, |Grimth, Yelle & Marley 


1998 


). The only published photometry for this object 



is that by Matthews et al. 1996, hereafter MNKO. 



Photometry is vital for determining Gliese 229B's en- 
ergy distribution (and hence luminosity and radius), 
as well as for flux calibrating the observed spectra. 

In this paper we present new JHKL' photometry 
for Gliese 229B which is more accurate than the pho- 
tometry obtained by MNKO. We also re-flux cali- 
brate the published spectra for the object and show 
that previous flux calibrations are in error. 

2. Observations 



We obtained JHKL' photometry for Gliese 229B 
on 1999 January 14 using the NASA Infrared Tele- 
scope Facility on Mauna Kea in Ha waii. We used the 
facility infrared camera NSFCAM flShure et al. 1994 ) 
together with the coronographic instrument CoCo 
dWang et al. 1994 [Toomey et al. 1998| ). CoCo is a 
cold coronographic attachment to NSFCAM that uti- 
lizes an apodized focal plane mask and a pupil mask 
that stops the telescope outer diameter down to 80% 
of the full aperture and also obscures the inner 30%. 
The focal plane mask is about 2.5 arcseconds in di- 
ameter. CoCo has an x/y adjustable pupil stop mask 
and a pupil imager for alignment. Flux in the wings 
of a stellar profile are reduced by a factor of 5 — 10 
when using CoCo. 

We also made use of the new tip-tilt capability of 
the telescope. The combination of the good image 
quality and tracking provided by the telescope with 
tip-tilt, and good removal of the light from the pri- 
mary star Gliese 229A by CoCo, made our images 
much cleaner than those obtained by MNKO, with a 
resulting increase in photometric accuracy. 



The observational technique used was to nod the 
telescope between Gliese 229A and a sky position 1 ar- 
cminute away. The individual exposure time on Gliese 
229B (and sky) at each of JHK was 90 seconds, and 
at L' it was 20 seconds. The total on-target integra- 
tion time was 6 minutes at J, 7.5 minutes at H, 7.5 
minutes at K and 20.7 minutes at L' . Two JHK flux 
standards were also observed (UKIRT faint standards 
4 and 12), and one V standard (HD 40335), with ap- 
propriate exposure times. The standards were placed 
on the same region of the detector as Gliese 229B, 
and were also observed by nodding to sky on alter- 
nate frames. 

The data were reduced by creating flatfields for ev- 
ery target and filter from the appropriate sky frames. 
Individual flats had to be created as the mask could 
flex slightly at different telescope positions and hence 
change the overall ilumination. At JHK each target 
frame was sky-subtracted and flatfielded. Aperture 
photometry was carried out on each frame and the er- 
ror taken to be the rms deviation in the set of Gliese 
229B observations, combined with the scatter in the 
standard star values. The combined error is 5% at 
JHK. At 1/ we obtained two sets of data, separated by 
a standard star observation. All the sky-subtracted 
and flatfielded frames in each set were combined (us- 
ing sigma-clipping in IRAF). The error was derived 
from the scatter in the remnant background levels, 
and is 10%. 

Our JHK values agree with those of MNKO within 
the errors, but our V value is three times brighter 
than their quoted value, which has very large errors. 
Our value agrees well with model predictions (see 
§4). Table 1 gives our results, and Figure 1 shows 
our reduced images at J (a single 90 second sky- 
subtracted and flatfielded image) and L' (the com- 
bined sky subtracted and flatfielded image for one of 
our two datasets, a total of 10 minutes on-target inte- 
gration). In Table 1 JHKL 7 are given on the UKIRT 
system, which is identical to the IRTF system except 
at H such that 

(H - K) UKIRT = 0.82[±0.02] x(H- K) IRTF 

flLeggett, Allard fc Hauschildt 1998b - 

3. Flux Calibrated Spectra and Luminosity 

We have recalibrated the published spectra for 
Gliese 229B using our new photometry and the pro- 
files of the NSFCAM filters. The filter profiles are for 
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Table 1 

PHOTOMETRY FOR GLIESE 229B 



J 


H K 


L' 


M-m* 


Luminosity 


BC K 




NSFCAM system 




mag 


logi O (L/L ) 


mag 


14.32 


14.35 14.42 


12.18 


1.185 


-5.18 


2.09 


±0.05 


±0.05 ±0.05 


±0.10 


±0.070 


±0.04 


±0.10 



van Altena, Lee, & Hoffleit 1994 







Q 


tra arose from their assumption that the target flux 
at the central wavelength of the filter could be simply 
scaled by the calibrator flux at that same wavelength. 
This is not a valid assumption as the flux distribu- 
tion across the filter is very different for Gliese 229B 
and any calibrator star. In the case of the Geballe 


T 






et al. 1996| spectra, those authors did not attempt 


m 






to scale by the observed photometry and their error 
arises from the usual problems of flux calibrating nar- 



Fig. 1.— NSFCAM and CoCo images of Gliese 229B 
at J (left) and L' (right). The image display stretch 
is zero to the Gliese 229B flux level. 



the operating temperature of the camera, and we have 
included the effect of the Mauna Kea atmospheric 
transmission. We have integrated the spectra over the 
JHKL' filters, and ratioed the flux to that obtained by 



integrating the observed Vega spectrum over the same 
filters. Vega was adopted to be zero magnitude at all 
wavelengths, and we scaled the Gliese 229B spectra 
such that the flux ratio through the filters matched 
our photometry. The atmosphere-convolved NSF- 
CAM filter profiles, and the Vega spectrum, are avail- 
able on request to S. Leggett (skl@jach.hawaii.edu), 
for calibration of other spectra. 

We found that although our JHK values are dif- 
ferent by <10% the published spectra in this region 
had to be multiplied by factors of 0.68 to 0.81, and al- 
though the L' magnitude differed by a factor of three, 
the published spectra in that region only had to be 
multiplied by a factor of 1.03. The error in the original 
flux calibration of the Oppenheimer et al. 1998 spec- 



row slit spectra combined with the uncertainty in re- 
moving the scattered light from Gliese 229A. 

Figure 2 shows the |Oppenheimer et al. 1998| spec- 
tra for Gliese 229B. We show the previously pub- 
lished spectral distribution, as well as our revised dis- 
tribution. The NSFCAM filter profiles are overlaid 
to demonstrate the unusual flux weightin g of the fil- 
ters for Gliese 229B. Figure 3 shows the |Geballe et 



al. 1996| higher resolution 1 — 2.5/mi spectrum, before 
and after recalibration. Scattered light from Gliese 
229A was seen in the raw data of both groups, despite 
their use of narrow slits. However the lack of flux in 
the strong methane- and water-absorption features 
shows that any residual light from Gliese 229A has 
been accurately removed. 

We have calculated the integrated luminosity us- 
ing the recalibrated spectra. The unobserved region 
between the K-band and L-band was represented by 
a linear interpolation, which is consistent with the 
models. The L-band spectrum ends at 4.15/im. The 
integrated flux beyond this point was adopted to be 
29% of the total integrated luminosity, as implied by a 
recent model by Allard & Hauschildt (private commu- 
nication). The luminosity in this region is affected by 
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Fig. 2. — Gliese 229B spectrum from pppenheimer 
et ill. 1998 — solid line shows the revised flux cal- 



ibration, dashed as published. The NSFCAM filter 
profiles are also shown. 




Wavelength yum 



Fig. 3. — Gliese 229B spectrum from Geballe et al 
19S|6| — solid line shows the revised flux calibration, 
dashed as published. 



strong CH4 and H2O bands at 6 — 8/im (unobservable 
from the ground), and this flux contribution is 25% 
smaller than that calculated by extending a Rayleigh- 
Jeans tail from the L spectrum. The total flux is de- 
termined to be 6.39e-15 W/m 2 with a 9% error (due 
to the flux calibration errors and an estimated 20% 
uncertainty in the flux beyond the L-band). This flux 
gives the luminosity and bolometric correction values 
listed in Table 1. MNKO use models to estimate a 
bolometric correction to their total in-band observed 
luminosity, and derive a total luminosity 3% smaller 
than ours (with an estimated 12% error). 

4. Comparison to Models 

Burrows et al. 1997 have published evolutionary 
models for giant planets and brown dwarfs. This work 



includes the condensation of grains in the equation of 
state but does not include opacity due to grains, and 
it is to be expected that these models will be revised 
shortly. Nevertheless our JHKL' colors and absolute 
magnitudes show good agreement with these models 
in the color-color diagrams (Burrows et al., Figures 
20 — 26) at all bands except the H-band. The discrep- 
ancy at H is large — about 0.4 magnitudes. However 
our three times brighter 1/ measurement (5 mag = 1.2) 
now gives a result consistent with the J and K values, 
implying that Gliese 229B has a mass between 25 Mj 
and 35 Mj and an age between 0.5 Gyr and 1.0 Gyr. 
The luminosity determined above and the Burrows et 
al. evolutionary tracks (their Figure 11) imply that 
a solution near the younger less-massive end of this 
range is more probable, and that the effective tem- 
perature is then around 900 K. (MNKO also deter- 
mined a temperature of 900 K through a comp arison 
of their pho tometry with a dust-free model by Tsuji| 



et al. 1996 



Burrows et al. suggest that the previous disagree- 
ment between the calculated and observed H and 1/ 
magnitudes were due to an incomplete CH4 opac- 
ity database. However CH4 is an important opac- 
ity source at K and L' both of which are now well 
matched by the models. Discrepancies are seen in the 
H-b and with the hotter M- and L-dwar fs also (see 
e.g. Leggett, Allard fc Hauschildt 1998 ) and these 
have been blamed on incomplete H2O opacity tables. 
Presumably the problem is due to incomplete knowl- 
edge of an opacity source, and the models, which are 
making large and rapid advances, will soon be able to 
match the entire Gliese 229B spectrum. 

5. Conclusions 

We have obtained new JHKL' photometry using 
the coronographic instrument C0C0 with NSFCAM 
on the IRTF. The data are more accurate than pre- 
viously published values, and while they agree within 
10% at JHK, we have determined an L' magnitude 
three times brighter than that previously published. 

We have re-flux-calibrated the published spectra 
for Gliese 229B by integrating the spectra over the 
NSFCAM filter profiles. We find that corrections of 
about 30% are required at JHK. The revised spectra 
are available in electronic form on request to B. Op- 
penheimer (bro@astro.caltech.edu) and T. Geballe 
(tgeballe@gemini.edu) . 

The JKL' colors and revised luminosity agree well 
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wit h the evolutio nary models of brown dwarfs by Bur- 
ro ws et aL 1997 , implying that Gliese 229B is likely 
to be a 0.5 Gyr-old 25 Mj object with T e// -900 K. 
We have not carried out a spectroscopic analysis as 
the models and synthetic spectra are currently being 
upgraded to include the effects of grains in both the 
equation of state and opacity tables — and grains 
have a large effect on the energy distribution of such 
cool objects. We can look forward to a better under- 
standing of the physics of atmospheres at these cool 
temperatures, now that there exists a more accurately 
flux calibrated spectrum for Gliese 229B. 

We are very grateful to the staff at IRTF. The 
IRTF, the NASA Infrared Telescope Facility, is op- 
erated by the University of Hawaii under contract to 
NASA. We are also grateful to France Allard and Pe- 
ter Hauschildt for calculating the flux contribution 
longwards of 4.15/im, and to the referee for useful 
comments. 
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